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Abstract: A sizeable level of non-Gaussianity in the primordial cosmological perturba- 
tions may be induced by a large trispectrum, i.e. by a large connected four-point correlation 
function. We compute the effect of a primordial non-Gaussian trispectrum on the halo mass 
function, within excursion set theory. We use the formalism that we have developed in a 
previous series of papers and which allows us to take into account the fact that, in the 
presence of non-Gaussianity, the stochastic evolution of the smoothed density field, as a 
function of the smoothing scale, is non-mar kovian. In the large mass limit, the leading- 
order term that we find agrees with the leading-order term of the results found in the 
literature using a more heuristic Press-Schecther (PS)-type approach. Our approach how- 
ever also allows us to evaluate consistently the sub leading terms, which depend not only 
on the four-point cumulant but also on derivatives of the four-point correlator, and which 
cannot be obtained within non-Gaussian extensions of PS theory. We perform explicitly 
the computation up to next-to-leading order. 
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1. Introduction 

Over the last decade a great deal of evidence has been accumulated from the Cosmic 
Microwave Background (CMB) anisotropy and Large Scale Structure (LSS) spectra that 
the observed structures originated from seed fluctuations generated during a primordial 
stage of inflation. While standard single-field models of slow-roll inflation predict that 
these fluctuations are very close to gaussian (see fl], |2|), non-standard scenarios allow 
for a larger level of non-Gaussianity (NG) (see Q and references therein). A signal is 
gaussian if the information it carries is completely encoded in the two-point correlation 
function, all higher connected correlators being zero. Deviations from Gaussianity are 
therefore encoded, e.g., in the connected three- and four-point correlation functions which 
are dubbed the bispectrum and the trispectrum, respectively. A phenomenological way 
of parametrizing the level of NG is to expand the fully non-linear primordial Bardeen 
gravitational potential $ in powers of the linear gravitational potential 3>l 

$ = $l + /nl ($l - (®b) + 3nl$nl • (1-1) 

The dimensionless quantities /nl and <?nl set the magnitude of the three- and four- 
point correlation functions, respectively [[|. If the process generating the primordial non- 
Gaussianity is local in space, the parameter /nl and <?nl in Fourier space are independent 
of the momenta entering the corresponding correlation functions; if instead the process 
which generates the primordial cosmological perturbations is non-local in space, like in 
models of inflation with non-canonical kinetic terms, /nl and <7nl acquire a dependence 
on the momenta. It is clear that detecting a significant amount of non-Gaussianity and its 
shape either from the CMB or from the LSS offers the possibility of opening a window into 
the dynamics of the universe during the very first stages of its evolution. Non-Gaussianities 
are particularly relevant in the high-mass end of the power spectrum of perturbations, i.e. 
on the scale of galaxy clusters, since the effect of NG fluctuations becomes especially visible 
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on the tail of the probability distribution. As a result, both the abundance and the clus- 
tering properties of very massive halos are sensitive probes of primordial non-Gaussianities 
|, 1, H 0, || i 0, 0, and could be detected or significantly constrained by the various 
planned large-scale galaxy surveys, both ground based (such as DES, PanSTARRS and 
LSST) and on satellite (such as EUCLID and ADEPT) see, e.g. @ and @. Further- 
more, the primordial non-Gaussianity alters the clustering of dark matter halos inducing 
a scale-dependent bias on large scales 14, 15, 16] while even for small primordial 



non-Gaussianity the evolution of perturbations on super-Hubble scales yields extra con- 
tributions on smaller scales [17, [Tq |. The strongest current limits on the strength of local 
non-Gaussianity set the /nl parameter to be in the range —4 < /nl < 80 at 95% confidence 
level 111 



While the literature on NG has vastly focussed on the impact on observables induced 
by a non-vanishing bispectrum, only recently attention has been devoted to the impact 



of a non- vanishing trispectrum of cosmological perturbations [2C, 21, 22, 23, 24]. This 
has been computed in several models, such as multifield slow-roll inflation model [25, 
26, 27, 28], the curvaton model p9f| , theories with non-canonical kinetic terms both in 
single field |3 
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3, 33] and in the multifield case [34], and in the case in which the 
cosmological perturbations are induced by vector perturbations of some non-abelian gauge 
field p5j| . While the most natural value of the <?nl parameter is 0(f^ L ), there are cases 
in which |#nl| 3> 1 even if /nl is tiny [36]. The effects of a cubic correction to the 
primordial gravitational potential onto the mass function and bias of DM haloes have been 



recently analyzed in [37] where the theoretical predictions have been compared to the results 
extracted from a series of large A— body simulations. The limit — 3.5-10 5 < <?nl < +8.2- 10 5 
has been obtained at 95% confidence level in the case in which the NG is of the local type. 

The goal of this paper is to present the computation of the DM halo mass function 
from the excursion set theory in the presence of a trispectrum, thus extending our previous 
computation of the DM halo mass function when the NG is induced by a bispectrum 
The halo mass function can be written as 



dn(M) 
dM 



p dlna^iM) 
~W d\nM 



X2) 



where n(M) is the number density of dark matter halos of mass M, a{M) is the variance 
of the linear density field smoothed on a scale R corresponding to a mass M, and p 
is the average density of the universe. Analytical derivations of the halo mass function 
are typically based on Press-Schechter (PS) theory |^9| and its extension [^, 41 1 known 



as excursion set theory (see [42] for a recent review). In excursion set theory the density 
perturbation evolves stochastically with the smoothing scale, and the problem of computing 
the probability of halo formation is mapped into the so-called first-passage time problem 
in the presence of a barrier. 

The computation of the effect of a primordial trispectrum on the mass function has 
been performed in || 43, 37], and is based on NG extensions of Press-Schechter theory 
|39j| . Performing the computation using the two different approximations proposed in || 
and |43]|, respectively, to evaluate the impact of a non- vanishing bispectrum on the DM 
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halo mass functions, one finds that, to leading order in the large mass limit, the predicted 
halo mass functions are the same for the two methods, but differ in the subleading terms, 
i.e. in the intermediate mass regime. Apart from understanding what is the correct result 
for the subleading terms, there is yet a fundamental reason why we wish to apply the 
excursion set theory to compute the DM mass function with a non-vanishing trispectrum. 
The derivation of the PS mass function in [[0]] requires that the density field S evolves with 
the smoothing scale R (or more precisely with the variance S(R) of the smoothed density 
field) in a markovian way. Only under this assumption one can derive the correct factor of 
two that Press and Schechter were forced to introduce by hand. As we have discussed at 



length in [44], this markovian assumption is broken by the use of a filter function different 
from a sharp filter in momentum space and, of course, it is further violated by the inclusion 
of non-Gaussian corrections. The non-markovianity induced by the NG introduces memory 
effects which have to be appropriately accounted for in the excursion set. As we will see, the 
mass function indeed gets "memory" corrections in the intermediate mass regime, which 
depend on derivatives of the correlators, and therefore cannot be computed with the NG 
extensions of PS theory studied in || 43 1. 



This paper is organized as follows. In Section £| we recall the basic points of the 
formalism developed in [44] for gaussian fluctuations, and extended in [38] to the NG case. 
In Section || we compute the NG corrections with the excursion set method induced by 
a large trispectrum, and we present our results for the halo mass function. Finally, in 
Section || we present our conclusions. 

2. The basic principles of the computation 
2.1 Excursion set theory 

In excursion set theory one considers the density field 5 smoothed over a radius R with 
a tophat filter in coordinate space, and studies its stochastic evolution as a function of 



the smoothing scale R. As it was found in the classical paper [41], when the density 5(R) 
is smoothed with a sharp filter in momentum space, and the density fluctuations have 
gaussian statistics, the smoothed density field satisfies the equation 

-QS-=<S), (2.1) 

where S = & 2 (R) is the variance of the linear density field smoothed on the scale R and 
computed with a sharp filter in momentum space, while f](S) is a stochastic variable that 
satisfies 

{ri(S 1 )r,(S 2 ))=6 D (S 1 -S 2 ) t (2.2) 



where 5d denotes the Dirac delta function. Equations ( |2.1| ) and (2.2) are the same as 
a Langevin equation with a Dirac-delta noise rj(S), with the variance S formally playing 
the role of time. Let us denote by H(<5, S)d5 the probability density that the variable 
S(S) reaches a value between 5 and 5 + d5 by time S. A textbook result in statistical 
physics is that, if a variable 5(S) satisfies a Langevin equation with a Dirac-delta noise, 
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the probability density 11(5, S) satisfies the Fokker-Planck (FP) equation 

dU 1 d 2 U 



(2.3) 



dS 2 d5 2 ' 

The solution of this equation over the whole real axis — oo < 5 < oo, with the boundary 
condition that it vanishes at 5 = ±oo, is 

n ° ( ^ ) = 7S e ~ 5V(25) - {2A) 

and is nothing but the distribution function of PS theory. Since large R, i.e. large halo 
masses, correspond to small values of the variance S, in ref. [41| it was realized that we 
are actually interested in the stochastic evolution of 5 against S only until the "trajectory" 
crosses for the first time the threshold 5 C for collapse. All the subsequent stochastic evo- 
lution of 5 as a function of S, which in general results in trajectories going multiple times 
above and below the threshold, is irrelevant, since it corresponds to smaller-scale structures 
that will be erased and engulfed by the collapse and virialization of the halo corresponding 
to the largest value of R, i.e. the smallest value of S, for which the threshold has been 
crossed. In other words, trajectories should be eliminated from further consideration once 
they have reached the threshold for the first time. In ref. jll]] this is implemented by 
imposing the boundary condition 

11(^)1^ = 0. (2.5) 
The solution of the FP equation with this boundary condition is 

U(S, S) = -L= \e~ s2 ^ - e -^ c -s?/(2S)-\ (2 6) 

and gives the distribution function of excursion set theory. The first term is the PS result, 
while the second term in eq. (|2.6| ) is an "image" gaussian centered in 5 = 28 c . Integrating 
this Il(<5, S) over d5 from — oo to 5 C gives the probability that a trajectory, at "time" S, has 
always been below the threshold. Increasing S this integral decreases because more and 
more trajectories cross the threshold for the first time, so the probability of first crossing 
the threshold between "time" 5 and 5 + dS is given by J r (S)dS, with 

HS) = -§g f C J5Ii{5-S). (2.7) 

With standard manipulations (see e.g. [42] or El]) one then finds that the function /(<t) 
which appears in eq. (L2) is given by 

f(a) = 2a 2 T(a 2 ), (2.8) 

where we wrote S = a 2 . Using eq. (|2.6|) one finds 



m = [-X 12 -*- 5ll ^\ (2.9) 



Observe that, when computing the first-crossing rate, the contribution of the gaussian 
centered in 5 = and of the image gaussian in eq. (|2.6| ) add up, giving the factor of two 
that was missed in the original PS theory. 
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2.2 Refinements of excursion set theory 



While excursion set theory is quite elegant, and gives a first analytic understanding of 
the halo mass function, it suffers of two important set of problems. First, it is based on 
the spherical (or ellipsoidal) collapse model, which is a significant oversimplification of 
the actual complex dynamics of halo formation. We have discussed these limitations in 
detail in [|4q| , where we also proposed that some of the physical complications inherent 
to a realistic description of halo formation can be included in the excursion set theory 
framework, at least at an effective level, by taking into account that the critical value for 
collapse is not a fixed constant S c , as in the spherical collapse model, nor a fixed function 
of the variance a of the smoothed density field, as in the ellipsoidal collapse model, but 
rather is itself a stochastic variable, whose scatter reflects a number of complicated aspects 
of the underlying dynamics. The simplest implementation of this idea consists in solving 
the first-passage time problem in the presence of a barrier that performs a random walk, 
with diffusion coefficient Db, around an average value given by the constant barrier of 
the spherical collapse model (more generally, one should consider a barrier that fluctuates 
over an average value given by the ellipsoidal collapse model). In this simple case, we 
found in |45|] that the exponential factor in the Press-Schechter mass function changes 
from exp{— 5 2 /2a 2 } to exp{— a5 2 /2a 2 }, where a = 1/(1 + Db)- In this approach all our 
ignorance on the details of halo formation is buried into Db- The numerical value of Db, 
and therefore the corresponding value of a, depends among other things also on the details 
of the algorithm used for identifying halos (e.g. on the link-length in a Friends-of-Friends 
algorithm). Observe that the replacement of exp{— 5 2 /2a 2 } with exp{— a5 2 /2a 2 } (with a 
taken however as a fitting parameter) is just the replacement that was made in refs. [^6|, , 
in order to fit the results of iV-body simulations. 

The second set of problems of excursion set theory is of a more technical nature, and 
is due to the fact that the Langevin equation with Dirac-delta noise, which is at the basis 
of the whole construction of ref. [41], can only be derived if one works with a sharp filter 
in momentum space, and if the fluctuations are gaussian. However, as it is well known, 
and as we have discussed at length in ref. [44], with such a filter it is not possible to 
associate a halo mass to the smoothing scale R. A unambigous relation between M and 
R is rather obtained with a sharp filter in coordinate space, in which case one simply has 
M = (4/3)7ri? 3 /9, where p is the density. When one uses a sharp filter in coordinate space, 
the evolution of the density with the smoothing scale however becomes non-markovian, 
and therefore the problem becomes technically much more difficult. In particular, the 
distribution function 11(5, S) no longer satisfies a local differential equation such as the FP 
equation. The issue is particularly relevant when one wants to include non-Gaussianities 
in the formalism, since again the inclusion of non-Gaussianities renders the dynamics non- 
markovian. 



In refs. [44, [3g] we have developed a formalism that allows us to generalize excursion set 
theory to the case of a non-markovian dynamics, either generated by the filter function or by 
primordial non-Gaussianities. The basic idea is the following. Rather than trying to derive 
a simple, local, differential equation for H(5, S) (which, as we have shown in ref. [44], is 
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impossible; in the non-markovian case II(<5, S) rather satisfies a very complicated equation 
which is non-local with respect to "time" S), we construct the probability distribution 
11(5, S) directly by summing over all paths that never exceeded the threshold S c , i.e. by 
writing II(<5, S) as a path integral with boundaries. To obtain such a representation, we 
consider an ensemble of trajectories all starting at So = from an initial position 8(0) = So 
and we follow them for a "time" S. We discretize the interval [0, S] in steps AS = e, so 
Sk = ke with k = l,...n, and S n = S. A trajectory is then defined by the collection 
of values {Si,...,S n }, such that 5(Sk) = 8k- The probability density in the space of 
trajectories is 

W(S ;S 1 ,...,5 n ;S n ) = (fe(£(Si) - S x ) ... 8 D (8(S n ) - 8 n )) , (2.10) 

where 8d denotes the Dirac delta. Then the probability of arriving in 5 n in a "time" S n , 
starting from an initial value So, without ever going above the threshold, is 1 

n e (J ; 5 n - S n ) = t dS x ... f° tW„_i W(S ; St,..., S n ^,S n ; S n ). (2.11) 

J —oo J —oo 

The label e in Il e reminds us that this quantity is defined with a finite spacing e, and 
we are finally interested in the continuum limit e — > 0. As we discussed in [^4|, |38|| , 
S\, ■ ■ ■ , Sn— i) Sn] S n ) can be expressed in terms of the connected correlators of the 
theory, 

W(S ;S 1 ,...,5 n ;S n ) = J V\ (2.12) 

{n oo / ,,p n n 

i + ^2 — — Yj---Yj Xi ^--- Xi p foi ■ ■ ■ K) 
i=l p=2 H=l ip=l 



where 



oo 



(2.13) 



Si = 5(Si), and (<5i . . . S n ) c denotes the connected n-point correlator. So 

n £ (<5 ; S n ; S n ) = f C dSt... dS n ^ I V\ (2.14) 

J —oo J 



exp < 



n 00 (—i)P n n 

ii=l p=2 P ' 



When S(S) satisfies eqs. (|2.1| ) and (2.2) (which is the case for sharp filter in momentum 



space) the two-point function can be easily computed, and is given by 

(S(S i )S(S j )) = mm(S i ,S j ). (2.15) 



In eqs. (2.4) and (2.6) we implicitly assumed So = 0. In the following however it will be necessary to 



keep track also of the initial position So. 
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If furthermore we consider gaussian fluctuations, all n-point connected correlators with 
n > 3 vanish, and the probability density W can be computed explicitly, 

W^(S ; 6 1 ,...,5 n ; S n ) = exp j-1 ]T(<5 i+1 - 5,) 2 } , (2.16) 



i=0 



where the superscript "gm" (gaussian-markovian) reminds that this value of W is com- 
puted for gaussian fluctuations, whose dynamics with respect to the smoothing scale is 
markovian. Using this result, in [44] we have shown that, in the continuum limit, the 
distribution function II e= o(<5; S), computed with a sharp filter in momentum space, satis- 
fies a Fokker-Planck equation with the boundary condition U e= o(5 c , S) = 0, and we have 
therefore recovered, from our path integral approach, the distribution function of excursion 
set theory, eq. Q2.6D - 

When the consider a different filter, eq. ( p. 15 ) is replaced by 

(SfaWSj)) = mm(S h Sj) + A(S h Sj) , (2.17) 

where A(Si,Sj) describes the deviations from a markovian dynamics. For instance, for a 
sharp filter in coordinate space, which is the most interesting case, the function A(Si,Sj) 
is very well approximated by 

A(Si, Sj) ~ k ~ — , (2.18) 

(for Si < Sj, and is symmetric under exchange of Si and Sj), with k ~ 0.45. The non- 
markovian corrections can then be computed expanding perturbatively in k. The com- 
putation, which is quite non-trivial from a technical point of view, has been discussed in 
great detail in [44]. Let us summarize here the crucial points. First of all, expanding to 
first order in A y - and using Aje*^ fcAfe<5ft = —idie l ^ k XkSk , where di = d/dSi, the first-order 
correction to IL is 



nf 1 (6 ; 5 n ; S n ) = dSi... d£ n _i - ^ A^didj 
J -°° ij=i 

/I n \ n 
VX exp < i ^i^i 5^ niin(5j, Sj)\i\j 

{ i=l i,3=l 

- &v / dSx... dSn-i didj W sm {S ; <5i , . . . , 5 n ; S n ) , 



(2.19) 



2 i,j=i - -» 

where we use the notation A^ = A(Si, Sj). One then observes that the derivatives di run 
over i = 1, . . . , n, while we integrate only over <i<5i . . . d5 n -\. Therefore, derivatives di with 
i = n can be simply carried outside the integrals. Derivatives di with i = 1, . . . , n — 1 are 
dealt as follows. Consider for instance the terms with i < n and j = n (together with 
j < n and i = n, which gives a factor of two). These are given by 

"- 1 r5 c 

J2 ^indn / d5i... dSn-i d l W sm (5 ; S h ..., S n ; S n ) , (2.20) 
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To compute this expression we integrate di by parts, 
f C dS x ... d5 n -t diW^ m (5 ; Si, ■ ■ ■ , $n, S n 

J — oo 



(2.21) 



dSi ...dSi. . .dS n -iW(S ;Si, . ..,Si = S c ,... ,S n -i,S n ;S n ) , 



where the notation dSi means that we must omit dSi from the list of integration variables. 
We next observe that W gm satisfies 



W gm (5 ; 5x,...,8i = 5 c ,... ,S n ; S n ) 

= W gm (S ; 5^, 8 C ; Si)W**(6 e ; 5 i+1 , ...,S n ;S n - Si) , 
as can be verified directly from its explicit expression Then 

d8\... dSi-i / d5 i+ i . . . d5 n -i 

-oo J — oo 

x^ m (,5o; Si-uSj Si)W^(S c ; S i+1 , . . . , S n] S n - Si) 

= nf m (5 ; S c ; S^ILf 1 ^; S n ; S n - Si) , 

and to compute the expression given in eq. ( 2.2C| ) we must compute 



(2.22) 



n-l 



£ A^nf m (5 ; S c ; Si)nf m (5 c ; S n ; S n - Si 



(2.23) 



(2.24) 



i=l 



To proceed further, we need to know Ilf m (<5o; S c ; Si). By definition, for e = this quantity 
vanishes, since its second argument is equal to the the threshold value S c , compare with 
eq. (|2.5| ). However, in the continuum limit the sum over i becomes 1/e times an integral 
over an intermediate time variable Si, 



n—l 



1=1 



S fi 



dSi , 



(2.25) 



so we need to know how Hf m (So; S c ; Si) approaches zero when e —* 0. In [44] we proved 
that it vanishes as i/e, and that 



Ur(5 ;5 c ;S) = Ve 



1 S c - So 

5 3 / 2 ' 



-(8 C -S f/(2S) 



Similarly, for S n < S c 



nr(5 c ;5 n ;S) = V~e-^ 6 ' 



n e -(S c -S n f/(2S) 



+ 0(e). 



+ 0(e) 



(2.26) 



(2.27) 



Therefore the two factor \J~i from eqs. ( 2.26| ) and ( |2.27 ) produce just an overall factor of 
e that compensate the factor 1/e in eq. ( 2.251 ), and we are left with a finite integral over 
dSi. Terms with two or more derivative, e.g. didj, or di,djdk acting on W, with all indices 
k mailer than n, can be computed similarly, and have been discussed in detail in [44]. 
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3. Contribution of the trispectrum to the halo mass function 



The effect of non-Gaussianities can be computed similarly, expanding perturbatively eq. ( |2.14| ) 
in terms of the higher-order correlators. In []38| we have examined the three-point correla- 
tor, i.e. the bispectrum. Here we compute the effect of the trispectrum. 

If in eq. (2.14) we only retain the four-point correlator, and we use the tophat filter in 
coordinate space, we have 

IL_(5 ; 6 n ; S n ) = d5 x . . . cZ<5 n _i T>\ 

J —CO J 

x exp I iXiSi - -[min(Si, Sj) + A, ;V \,A ; + — (5i5j5 k 5i}\i\j\k\i > ■ (3.1) 
Expanding to first order, Ajj and {5i5j5kSi} do not mix, so we must compute 

n[ 4) (<5 ; 5 n ; S n ) = i- £ (5. Wz> t db x ... ^ n _i^d fe cW gm , (3.2) 

i,j,k,l=l 00 

where the superscript (4) in IL^ refers to the fact that this is the contribution linear in 
the four-point correlator. 

In principle this expression can be computed with the same technique discussed above, 
separating the various contributions to the sum according to whether an index is equal or 
smaller than n. In this way, however, the computations faces some technical difficulties. 2 
Fortunately, the problem simplifies considerably in the limit of large halo masses, which is 
just the physically interesting limit. Large masses mean small values of S n . In eq. ( |3.2j ) 
the arguments SijSjjSk and Si in the correlator (5i5jdk5i) = (5(Si)5(Sj)8(Sp.)5(Si)) range 
over the interval [0, S n ] and, if S n goes to zero, we can expand the correlator in a multiple 
Taylor series around the point Si = Sj = Sk = S[ = S n . We introduce the notation 



dP cP d r d s 
dSfdS q j dS r k dSf^"" ; 



(SiSiMSjWSkMSi)) 



(3.3) 



Si — Sj — S]$ — Si — Sn 



Then 



(#(W)«= f; ( (3 - 4) 

p,q,r,s=0 

x(s n - s t ns n - SjYiSn - s k ) r (s n - StyG$«> r ''\s n ) . 

Terms with more and more derivatives give contributions to the function f(cr), defined 
in eq. (|1.2|) , that are subleading in the limit of small a, i.e. for a/S c -C 1, and the 
leading contribution to the halo mass function is given by the term in eq. ( |3.4j) with 
p = q = r = s = 0. At next-to-leading order we must also include the contribution of the 



2 In particular, when we have terms with three or more derivatives, we need to generalize eq. (2.27), 
including terms up to C(e 3//2 ), which is quite non-trivial. 
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terms in eq. (|3.4j) with p + q + r + s = l, i.e. the four terms (p = 1, q = 0, r = 0, s = 0), 
(p = 0, g = 1, r = 0, s = 0), (p = 0, q = 0, r = 1, s = 0) and (p = 0, q = 0,r = 0, s = 1); 
at next-to-next-to-leading order we must include the contribution of the terms in eq. ( fO| ) 
with p + g + r = 2, and so on. For the purpose of organizing the expansion in leading term, 
subleading terms, etc., we can reasonably expect that, for small S n 

G { t q ^ s \s n ) ~ S-^ +r+s \5 4 (S n )) , (3.5) 

i.e. each derivative d/dSi, when evaluated in Si = S n , gives a factor of order 1/S n . This 
ordering will be assumed when we present our final result for the halo mass function below. 
However, our formalism allows us to compute each contribution separately, so our results 
below can be easily generalized in order to cope with a different hierarchy between the 
various G^'' 9 ' r (S n ). 

The leading term in is 

/J}4\ n rS c 

U^ L \5 ;5 n ;S n ) = ^- £ / dS 1 . . . dSn^d^W^ , (3.6) 

where the superscript "L" stands for "leading". Since in the end we are interested in the 
integral over d5 n of H e (5o; 5 n ; S n ), see eq. (2.7), we can write directly 

f C d5 n (OS S n] S n ) = {6i( f n)) r ■ ■ ■ dSndidjO^W^ . (3.7) 

J -°° i,j,k,l=l J ~°° 

This expression can be computed very easily by making use of identities that we proved in 
]4~i| . [38[ |. Namely, we consider the derivative of Tlf m with respect to the threshold S c (which, 
when we use the notation nf m ((5o; 5 n ;S n ), is not written explicitly in the list of variable on 
which nf m depends, but of course enters as upper integration limit in eq. ( p.ll| )). Then 
one can show that 

f dS x ... d5 n cW gm = J- f C d5 n nf m , (3.8) 

■ i J — OO Ode J—oo 



i=l 



d8 x ... d5 n didj I4^ m = — / d5 n nf m , (3.9) 

■ i */ — OO C " — OO 



and similarly for all higher-order derivatives, so in particular 



I " d5x . . . dS n didjd^W^ 1 = —j I C d5 n IFf 1 . (3.10) 

i,j,h,l=l °° ad cJ-oo 

Therefore, in the continuum limit, the right-hand side of eq. ( p.Tj ) is computed very simply, 
just by inserting in eq. JOcD the value of nf m for e -»■ 0, 

n^ (5 = 0; 5 n - S n ) = -±= \e</W - e -(24=-W(2S«)l > {3U) 



- 10 - 



and therefore, in the continuum limit, 

(Sc 



5 C [3 



S n 



(3.12) 



We now insert this result into eqs. ( |2.7D and ( |2.8| ) and we express the result in terms of the 
normalized kurtosis 

(3.13) 



1 



<> 4 (a)^<<5 4 (S)>. 

Putting the contribution of I[( 4 ' L ) together with the gaussian contribution, and writing 
S = a 2 , we find 



OX 1 / 2 A 



7T 



a 



1 + 



a 2 S i (a) 
24 



45 2 



-4-^-3 +-a 



1 



24' diner \a 2 



(3.14) 

Let us emphasize that the variance a 2 is to be computed applying the linear transfer 
function to the primordial gravitational potential (1.1) containing the extra g^L-piece. The 
result given in eq. (3.14) agrees with the one obtained in |43| by performing the Edgeworth 
expansion of a non-Gaussian generalization of Press-Schechter theory. However, just as we 



have discussed in []38| for the case of the bispectrum, eq. (3.14) cannot be taken as the full 
result beyond leading order. If we want to compute consistently to NL order, we need to 
include the terms with p + q + r + s = lin eq. (3.4), which is given by 

/ dS n n( 4 < NL )(5 ; 5 n ; S n ) = -± Gf'°'°'°\s n ) £(S n - Si) 



j,k,l=l 



d5 1 ... dSn-xdSndidjdkdiWB* , (3.15) 



where the superscript "NL" in ]le 4,NL ^ stands for next-to-leading, and we used the fact that 
the four terms (p = 1, q = 0,r = 0, s = 0), . . . , (p = 0, q = 0, r = 0, s = 1) give the same 
contribution. We now use the same trick as before to eliminate ^™ k l=1 djdkdk in favor of 
d 3 /d5 3 c , 

d5 n U^ L \5o; 5 n ; S n ) = -— G^ 0fi ' 0) (S n ) £(S„-Si)^ / d5 1 . . . d5 n ^d5 n diW^ . 



i=l 



(3.16) 

The remaining path integral can be computed using the technique discussed in eqs. fl2.21|) - 
Q2.27[), and we get 



d5 n U^(5 ; 5 n ; = G^'°' 0,0 \s n 

Z47T 



dSi 



sf /2 (S n - Si) 1 / 2 



S 2 J(2S i ) 



d5 n (5 C - 6 n ) exp 



2(5' J1 Si) 



(3.17) 



Observe that this expression involves an integral over all values of an intermediate "time" 
variable Si which ranges over < Si < S n . As we have discussed in detail in [44, 38] these 
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terms are "memory" terms that depend on the whole past history of the trajectory, and 
reflect the non-markovian nature of the stochastic process. 
The integral over d5 n is easily performed writing 

so it just gives (S n — Si). Carrying out the third derivative with respect to 5 C and the 
remaining elementary integral over dSi we get 



r s c * r f(i,o I o,o) r<? s 

/ dS n n(W(* ; 5n - S n ) = ^ J Sn) e~ s ^ . (3.19) 



We now define 

tU{*) = 4 s2 { } , (3-20) 



where as usual S = a 2 . When the ordering given in eq. (3.5) holds, U±(g) is of the same 



order as the normalized kurtosis 54(0"). Computing the contribution to f(cr) from eq. ( 3.1S| ) 
and we finally find 



1 , o*Si(v) (5t A5 2 \ 1 
i + ^4~~ V ~^ J 24 



24 \a 2 J 24 dlna 



d In a \a 2 

(3.21) 



The above result only holds up to NL order. If one wants to use it up to NNL order, 
the terms in square bracket where ^S^a), a 2 dS<i/dhia, o 2 U 4 (o), and a 2 dM<i/d In a are 
multiplied by factors of 0(1) must be supplemented by the computation of the terms 
with p + g + r + s = 2in eq. ( |3,4j ), which will give a contribution analogous to the term 
Vz{(j) computed in [ ]38| | . However, with present numerical accuracy, NNL terms are not yet 
relevant for the comparison with A^-body simulations with non-Gaussian initial conditions. 

Our result may be refined in several ways. Until now we have worked with a barrier 
with a fixed height 5 C and we neglected the corrections due to the filter. We can now 
include the modifications due to the fact that the height of the barrier may be thought 
to diffuse stochastically, as discussed in |4q| , and also the corrections due to the filter. 
To compute the non-Gaussian term proportional to the four-point correlator with the 



diffusing barrier we recall, from [45|, that the first-passage time problem of a particle 
obeying a diffusion equation with diffusion coefficient D = 1, in the presence of a barrier 
that moves stochastically with diffusion coefficient Db, can be mapped into the first-passage 
time problem of a particle with effective diffusion coefficient (1 + Db), and fixed barrier. 
This can be reabsorbed into a rescaling of the "time" variable S — > (1 + Db)S = S/a, 
and therefore a — > a/y/a. At the same time the four-point correlator must be rescaled 
according to (5%) — > a^ 2 (5^) since, dimensionally, (5%) is the same as S 2 , which means 
that ^4 — * aS 4 . As a final ingredient, we must add the effect of the tophat filter function 
in coordinate space. For a tophat filter in coordinate space, we have found [44] that the 
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two-point correlator is given by eqs. ( 2,17] ) and ( 2.18; ), Including the non-markovianity 



induced by the tophat smoothing function in real space, using the computations already 
performed in [44, 38 1, we end up with 

/(,) = (l-,)0) 1/2 ^ e -^) 

a 2 S 4 (a) (aHi _ 4a5 2 _ \ 1 2 dS 4 ( a% _ 
i+ 24 U 4 a 2 V 24 a dlna [ a 2 6 



a 2 U 4 (a) faSi \ 1 2 



24 U 2 / 24 dlncr 



+ 1 - ttt^ 



K a 



X / 2 <L / a5\ 



\2 



+ (3 - 22) 



where ft = k/{\ + Dg). This is our final result. More generally, also the term pro- 
portional to the incomplete Gamma function could get non-Gaussian corrections, which 
in principle can be computed evaluating perturbatively a "mixed" term proportional to 
^ij{bk$ibm,$n)didjdkdid m d n . However we saw in that in the large mass limit, where 
the non-Gaussianities are important, the term proportional to the incomplete Gamma 
function is subleading, so we will neglect the non-Gaussian corrections to this subleading 
term. 

4. Conclusions 

In this paper we have computed the DM halo mass function as predicted within the ex- 
cursion set theory when a NG is present under the form of a trispectrum. We thus have 
extended our previous results presented in ref. [38], where a similar computation was per- 



formed in the presence of a NG bispectrum. Our computation accounts for the non- 
markovianity of the random walk of the smoothed density contrast, which inevitably arise 
when deviations from gaussianity are present. While our result coincides at the leading 



order 0(5^/a 4 ) with that obtained in [43] through PS theory, it is different at the order 
0{5 2 c /a 2 ). This is due to the memory effects induced by the non-markovian excursion set 



which are not present in the PS approach. Our final expression ( 3.22j ) takes into account 



as well the non-markovian effects due to the choice of the tophat filter in real space and 
the proper exponetial decay at large DM masses. 

The work of MM is supported by the Fond National Suisse. The work of AR is supported 
by the European Community's Research Training Networks under contract MRTN-CT- 
2006-035505. 
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